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This  is  the  first  in  a series  of  reports  dealing  with  the  results  of  a major  DNA  program 
coordinating  rocket,  satellite,  and  ground-based  equatorial  studies  at  Kwajaiein  Atoll  in  the 
Marshall  Islands.  This  report  concentrates  on  the  NRL  contribution  to  the  overall  program,  the 
direct  measurement  of  the  ionospheric  plasma  state  as  determined  by  on-board  pulsed-plasma- 
probe  (P^)  observations  of  election  energies,  densities,  and  density  fluctuation  power  spectra. 

In  particular,  we  describe  our  newly  developed  data  acquisition  system  which  acquired  and 
displayed  in  real-time  the  ionospheric  electron  density  profile  as  it  was  being  measured  by  the 
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REAL-TIME  DISPLAY  OF  IONOSPHERIC  ELECTRON  DENSITY  PROFILES 
FROM  A ROCKET-BORNE  PLASMA  PROBE: 

The  Payload,  Acquisition  System,  and  the  Real-Time  Results 
I.  INTRODUCTION 

On  23  August  1977,  an  ionospheric  rocket  payload  was 
launched  from  Roi-Namur  in  the  Kwajalein  Atoll  to  study 
the  cause-effect  relationships  between  equatorial  plasma 
instabilities,  ionospheric  irregularities,  and  scintillation 
phenomena . 

The  payload,  illustrated  in  Figure  1,  was  in:  :rumented 
with  an  ion  mass  spectrometer,  two  photometers  (6300  k 
and  5577  k) , electric  field  sensors,  and  a number  of  probes 
which  utilized  various  techniques  for  the  determination 
of  electron  density,  temperature,  and  density  fluctuation 
power  spectra. 

The  launch,  as  part  of  a major  DNA  program  carrying  the 
code  name  EQUATORIAL  WIDEBAND,  was  coordinated  with  Altair 
radar,  top-  and  bottom-side  soundings,  and  ground-based 
photometric  measurements  of  F-region  winds.  These  correlative 
observations  permitted  near-continuous  monitoring  of  the 
ionospheric  plasma  state  with  emphasis  on  the  develop- 
ment and  morphology  of  regions  of  plasma  irregularities. 

The  timing  of  the  launch  was  critically  hinged  to 
space-time  coincidence  with  the  Wideband  Satellite  transit 


Note:  Manuscript  submitted  November  3,  1977. 


so  that  optimum  correlation  could  be  achieved  between  in- 
situ  observations  of  ionospheric  irregularities  and 
scintillation  spectra  as  measured  by  Wideband. 

This  report  will  focus  on  the  NRL  contribution  to 

the  overall  program,  the  direct  measurement  of  the 

ionospheric  plasma  state  as  determined  by  on-board  pulsed- 
3 

plasma-probe  (P  ) observations  of  electron  energies, 

densities,  and  density  fluctuation  power  spectra.  In 

particular  we  describe  our  newly  developed  data  acquisition 

system  which  acquired  and  displayed  in  real-time  the 

ionospheric  electron  density  profile  as  it  was  being 

3 

measured  by  the  P technique  during  the  rocket  flight. 

The  development  of  the  real-time  system  was  part 

of  the  NRL  effort  to  expedite  the  transition  from  data 

collection  to  information  transfer  and  publication.  The  l 

capability  was  also  developed  for  future  programs  involving 

ionospheric  modification  experiments,  where  in  some  cases 

the  placement  of  upper  atmospheric  chemical  releases 

can  be  too  critical  to  rely  on  predicted  trajectory 

information  and  ground-based  lonosonde  determinations 

of  the  F-layer  peak.  This  subject  is  discussed  in  the 

closing  section  of  this  report.  Before  that  we  describe 
3 

the  P experiment,  the  acquisition  and  display  system, 
and  the  real-time  results. 
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XI.  THE  NRL  PULSED  PLASMA  PROBE  (P3)  EXPERIMENT 

3 

Since  the  P technique  has  been  described  in  detail 
1-4 

elsewhere  , only  a brief  account  will  be  given  here 
in  order  to  facilitate  an  understanding  of  the  data  stream, 
acquisition  system,  and  the  relevance  to  geophysical  units. 

The  pulsed  plasma  probe,  shown  in  the  payload  sketch 
of  Figure  1,  carried  the  time-dependent  voltage  illustrated 
in  Figure  2.  A continuous  chain  of  sawtooth  sweeps  was 
automatically  applied  to  the  probe  and  its  guard  electrode, 
with  every  9th  and  10th  sawtooth  replaced  by  a constant 
dc  level  V_.  The  currents  I sampled  during  the  voltage 

o S 

pulses  V on  the  sawtooth  envelope  were  used  to  generate 
s 

conventional  current-voltage  ( I , V ) characteristics  from 

s s 

which  absolute  electron  density  and  temperature  are  determined 
The  interpulse  sampling  of  currents  (3  Ig)  at  the  fixed 
baseline  voltage  VD  tracked  ionospheric  density  variations 

O 

which  may  have  occurred  during  the  sweep  time  T , and 

S 

provided  the  raw  data  for  the  determination  of  electron 

density  fluctuation  power  spectra.  The  IQ  values  were 

sampled  at  a 512  and  1024  Hz  rate  during  the  pulsed  sweep 

and  fixed  bias  periods,  respectively. 

Nominal  values  for  sweep  voltage  parameters  (as 

defined  in  Fig.  2)  were: 

<V_,  V+)  - (-1.55v,  +3 . Ov) 

VB  - +2.0  v 

r - 400  ms 

s 


3 


The  electrometer  response  covered  5 full  decades 
of  current  with  a system  of  automatic  switching  which 
utilized  6 linear  electrometer  ranges.  Twelve  monitor 
levels  uniquely  identify  the  range  of  operation  and  the 
polarity  of  current  collection.  In  this  system  every  0 to  5 
volt  telemetry  value  for  the  electrometer  output  ^TMV(Ig) 
or  TMV(Ib)j  had  a range  value  (Rg  or  Rg)  which  determined 
the  absolute  magnitude  and  polarity  of  current  collection 
through  the  linear  relationship 

Vb  • an*  ™V<IS,B>  ♦ bn  <l) 

n - 1 12 

With  this  system,  "on  scale"  operation  is  always  achieved 

I 

for  every  I-V  characteristic,  regardless  of  the  ionospheric 
plasma  state. 

III.  THE  ACQUISITION  AND  DISPLAY  SYSTEM 
The  specific  tasks  accomplished  by  the  NRL  real- 
time data  acquisition  system  were: 

i)  The  acquisition  and  storage  on  disk  file  via  the 
PDP  11/10  digital  computer  all  data  relevant  to 
the  NRL  Pulsed  Plasma  Probe  (P“*)  in  the  PCM 
Telemetry  Stream  (i.e.,  Ig,  Ig , Vg,  Rg,  Rg  and 
magnetometer  response),  and 

4 

ii)  The  reduction  and  display  in  real-time  of  the 

relative  electron  density  profile  (indicated  by 


4 


baseline  current  measurements  IQ)  as  a function 
of  time  and  predicted  altitude. 

We  describe  below  in  general  terms  the  programming 
effort  involved  and  discuss  some  of  the  more  important 
details  in  a real-time  acquisition  system. 

The  general  method  employed  in  the  data  acquisition 
and  display,  is  illustrated  in  the  flowcharts  of  Figs. 

3-5.  There  are  two  Interrupt  Service  Routines  (ISR's) 
and  one  main  Fortran  routine.  The  sequence  is  driven 
by  interrupts  to  the  PDP  11/10  Digital  Data  Interface 
Unit  (DR11-K)  control  lines.  The  interrupts,  which  occur 
at  the  beginning  of  pre-selected  words,  were  supplied  by 
a CONIC  decommutator  unit  (DPAD) . The  unit  is  programmable 
and  facilitates  interrupt  position  changes  throughout  the 
full  32  word  range  of  a subframe*. 

A.  Raw  Data  Description 

Words  4,  6,  26,  and  27  in  the  PCM  format  shown  in 

3 

Fig.  6 were  relevant  to  the  P experiment.  Below  we 
describe  the  significance  of  each  of  the  main  frame  words. 

* 

NOTE:  One  mainframe  consisted  of  16  minor-frames  (or 
sub-frames)  each  of  32  words  (0-31)  in  length.  The  data 
rate  of  the  PCM  stream  was  1024  subframes/sec  - 64  main- 
frames/sec - 64  x 16  x 32  (32.768  x 10^)  words/sec.  Fig. 


illustrates  one  mainframe  and  the  position  of  relevant  NRL 
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P baseline  current  Ig  and  sweep  current  Ig  appeared 
alternately  in  word  26  of  each  subframe  (i.e.,  they  were 
subcommutated  at  a 512  Hz  rate).  The  associated  range 
monitors,  R„  and  R_  , followed  immediately  in  word  27. 

D S " 

In  addition,  in  every  8th  occurrence  of  word  27,  the 
Ig  range  monitor  value  was  replaced  by  the  applied  probe 
voltage  sweep  monitor  V . 

Word  6 in  subframes  3,  7,  and  11  of  each  mainframe 
provided  a ground  reference,  a 4-5  volt  reference,  and 
magnetometer  data,  respectively.  The  frequency  of  magnetom- 
eter data  sampling  was  therefore  once  per  mainframe 
(-  16  subframes)  or  64  Hz. 

Word  4 in  subframe  1 of  each  mainframe  was  used 
as  a synchronization  Indicator  for  the  TM  signal  and 
served  as  a flag  marking  lift-off. 

B.  Interrupt  Service  Routines 

DRROU  sorted  data,  and  set  up  and  filled  three  data 
buffers  using  a double  buffering  scheme.  Each  data  type 
(i.e.,  In,  I ,...)  was  uniquely  tagged  for  later  reference 
by  the  reduction  software.  The  three  buffers  handled  by 
this  routine  were  those  containing:  (1)  I0,  Ig , vg^g£p» 

(2)  Ig  (for  display  only),  and  (3)  magnetometer,  ground  and 
5 volt  reference. 

DRROU  performed  no  Input  or  Output  operation  but  set 
flags  for  either  ARROU  or  the  main  FORTRAN  routine  to  test 
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for  transfer  to  disk  or  display.  A buffer  was  full  when 
the  word  count  reached  512. 

ARROU  kept  time  (the  real-time  clock  was  programmed 
to  interrupt  on  overflow  every  millisec),  transferred  full 
data  buffers  to  disk,  and  tested  flags  set  by  DRROU  to 
determine  whether  the  ID,  I , or  magnetometer  buffers 
were  full. 

C.  Main  FORTRAN  Routine  (IBDIS3) 

IBDIS3  provided  for  initialization  and  plotting 
routine  set-up,  as  well  as  the  Ig  calculation  from  Eq  • (1). 
The  Ig  values  were  then  averaged  over  1.5  sec  intervals 
for  display  on  a Tektronix  video  controller. 

IV.  REAL-TIME  RESULTS 

In  the  application  reported  here,  consecutive  1.5  sec 
time  elements  were  blocked  in  real-time  and  the  averaged 
value  of  Ig  plotted  as  a measure  for  relative  electron 
density  variations.  Averaging  conserved  program  running 
time  and  smoothed  payload  spin  modulation  of  collected 
currents.  A replay  of  the  real-time  results  is  presented 
in  Figure  7 with  the  bottom  abscissa  showing  time  after 
lift-off  and  the  upper  abscissa  showing  the  actual  altitude. 
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Detailed  post-flight  determinations  of  electron 
density  through  standard  analysis  of  the  Ig-V,, 

characteristics  proved  the  real-time  <IB>  profile  to 

be  an  excellent  monitor  of  relative  electron  density, 

the  proportionality  being  <IQ>-  lO-4  A “ » 3 (104)cra~3. 

Reference  6 discusses  some  limitations  on  the  assumption 

I « ji 
1 B e 

We  note  that  the  payload  had  a predicted  apogee 
of  850  km  but  only  reached  an  altitude  of  160  km  as 
a result  of  a failure  in  second  stage  separation. 

The  fact  that  the  payload  never  got  to  the  F-layer  was 
evident  in  the  real-time  relative  density  profiles 
that  did  not  show  the  level  of  currents  expectea 

for  F-region  penetration. 

V.  COMMENTS  AND  CONCLUSIONS 

We  have  described  our  recently  developed  real- 
time data  acquisition  and  display  system  with  emphasis 
in  its  first  application  in  a rocket  program  funded 
by  the  Defense  Nuclear  Agency. 

In  real-time  the  system  acquired  the  PCM 
elements  of  a rocket  borne  pulse  plasma  probe  and  plotted 
the  relative  electron  density  profile  as  it  was  being 
measured  during  the  flight. 
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An  unfortunate  failure  in  the  vehicle's  second 
stage  separation  served  to  demonstrate  one  of  the  important 
features  of  the  real-time  system  particularly  as  it  can 
be  applied  to  future  programs  involving  chemical 

releases  at  ionospheric  altitudes.  The  vehicle  was  to 
deliver  the  payload  to  the  top-side  of  the  F-region 
but  only  achieved  a 160  km  apogee... a result  which  was 
evident  in  the  real-time  density  profile.  Had  an 
F-region  chemical  release  been  an  objective,  the  mission 
would  have  failed.  However,  future  programs  with  the 
capability  for  real-time  in-situ  profile  determination 
at  the  position  of  the  explosive  payload  could  easily 
pinpoint  primary  target  areas  (such  as  the  bottom-side 
gradient  of  the  F-layer  or  the  F-layer  peak)  as  well 
as  secondary  targets  (like  sporadic-E) . In  the  vehicle 
failure  described  here,  easily  considered  catastrophic 
if  an  F-region  release  was  intended,  a secondary  success 
coull  have  been  achieved  by  retining  the  release  for 
the  downleg  penetration  of  the  E-region  irregularities 
shown  in  Fig.  7. 

Some  Software  considerations.  There  were  two  areas 
which  required  special  attention  in  the  real-tine  da-a 
acquisition  and  display  routines.  Briefly,  they  were  in 
the  areas  of : 

i)  Synchronization  loss  (or,  signal/interrupt  loss)  and, 
ii)  Loop  timing. 
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These  two  elements  were  not  independent  in  that  the 
additional  software  necessary  to  account  for  signal 
losses  exceeded  the  allowed  loop  time  of  the  ISR  DRROU. 


Since  I0  and  Ig  occurred  alternately  in  word  26, 
the  acquisition  routines  had  to  track  and  identify  each 
subframe  as  it  was  sampled.  It  was  necessary  to  detect 
a loss  or  distortion  of  this  information  so  that  program 
logical  flows  were  not  corrupted.  We  note  that  in  addition 
to  the  interrupt  an  even  or  odd  subframe  indicator  would 
be  desirable.  This  was  not  originally  available  but  was 
finally  implemented  by  synchronization  through  interrupt 
counting . 

The  final  real-time  display  routine  was  run 
independently  of  the  permanent  data  acquisition  soft- 
ware Storage  to  the  PDP/11  disk  medium  was  accomplished 
immediately  after  the  flight  by  replaying  the  flight 
PCM  stream  from  the  flight  tape  at  one-half  the  actual 
data  rate  (16.38  7/sec).  Half  speed  playback  was 
necessitated  by  the  increased  loop  time  resulting  from 
synchronization  problems.  In  the  future  this  dual 
procedure  will  not  be  necessary  since  a simpler  method 
of  re-establishing  sync  will  be  implemented  through 
hardware.  Following  playback  and  within  one  hour  after 
flight  termination  , detailed  data  analysis  routines  (Nfi) , 

6N  , and  T ) were  initiated.  These  results  will  be  the 
e e 

subject  of  a follow-up  report. 
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Fig.  1 — DNA  equatorial  wideband  rocket  payload 
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Fig.  5 — DRROU  — DPAD  interrupt  service  routine 
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